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Abstract. In this article, we consider a generalized Radon trans- 
form that comes up in ultrasound reflection tomography. In our 
model, the ultrasound emitter and receiver move at a constant dis- 
tance apart along a circle. We analyze the microlocal properties of 
the transform 1Z that arises from this model. As a consequence, we 
show that for distributions with support sufficiently inside the cir- 
cle, 71*71 is an elliptic pseudodifferential operator of order —I and 
hence all the singularities of such distributions can be recovered. 



1. Introduction 

Ultrasound reflection tomography (URT) is one of the safest and 
most cost effective modern medical imaging modalities (e.g. see [9, 10, 
11, 12] and the references there). During its scanning process, acoustic 
waves emitted from a source reflect from inhomogeneities inside the 
body, and their echoes are measured by a receiver. This measured data 
is then used to recover the unknown ultrasonic reflectivity function, 
which is used to generate cross- sectional images of the body. 

In a typical setup of ultrasound tomography, the emitter and re- 
ceiver are combined into one device (transducer). The transducer emits 
a short acoustic pulse into the medium, and then switches to receiv- 
ing mode, recording echoes as a function of time. Assuming that the 
medium is weakly reflecting (i.e. neglecting multiple reflections), and 
that the speed of sound propagation c is constant 1 , the echoes measured 
at time t uniquely determine the integrals of the reflectivity function 
over concentric spheres centered at the transducer location and radii 
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This assumption is reasonable in ultrasound mammography, since the speed of 
sound is almost constant in soft tissue. 
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Figure 1 . A sketch of integrating curves in URT 



r = ct/2 (see Fig. 1 (a) below, [12] and the references there). By focus- 
ing the transducer one can consider echoes coming only from a certain 
plane, hence measuring the integrals of the reflectivity function in that 
plane along circles centered at the transducer location [11]. Moving 
the transducer along a curve on the edge of the body, and repeating 
the measurements one obtains a two-dimensional family of integrals of 
the unknown function along circles. Hence the problem of image re- 
construction in URT can be mathematically reduced to the problem 
of inverting a circular Radon transform, which integrates an unknown 
function of two variables along a two-dimensional family of circles. 

In the case when the emitter and receiver are separated, the echoes 
recorded by a transducer correspond to the integrals of the reflectivity 
function along confocal ellipses with foci corresponding to the locations 
of the emitter and receiver (see Fig. 1(b)). While this more general 
setup has been gaining popularity in recent years (e.g. see [9, 10]), 
the mathematical theory related to elliptical Radon transforms is rel- 
atively undeveloped. In this paper we study the microlocal properties 
of this transform which integrates an unknown function along a family 
of ellipses. 

2. Definitions and Preliminaries 

We will first define the elliptical Radon transform we consider, pro- 
vide the general framework for the microlocal analysis of this transform, 
and show that our transform fits within this framework. 

2.1. The Elliptical Transform. In URT the locations of the emitter 
and receiver are limited to a curve surrounding the support of the 
function to be recovered and the data taken can be modeled as integrals 
of the reflectivity function of the object over ellipses with foci being the 
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Figure 2. A sketch of the domain and the notations 

transmitter and the receiver. In this paper we consider this curve to be 
a circle, which is both the simplest case mathematically and the one 
most often used in practice. By using a dilation, we can assume the 
circle has radius r = 1. We also make the restriction that the source 
and detector rotate around the circle a fixed distance apart. We denote 
the fixed difference between the polar angles of emitter and receiver by 
2a, where a £ (0, 7r/2) (see Fig. 2) and define 

a = sina, b = cosa. (2-1) 

As we will see later our main result relies on the assumption that 
the support of the function is small enough. More precisely, we will 
assume our function is supported in the ball 

D b = {x e R 2 \\x\ < b}. 

We parameterize the trajectories of the emitter and receiver, respec- 
tively, as 

7t(s) — (cos(s — a), sin(s — a)) 

1r{s) = (cos(s + a), sin(s + a)) for s £ [0, 2%]. 

Thus, the emitter and receiver rotate around the unit circle and are 
always 2a units apart. For s £ [0, 2n] and L > 2a, let 

E(s, L) = {x £ R 2 | \x - 7 t0)| + \x- 7*0)1 = L}. 

Note that the center of the ellipse E(s,L) is (b cos s, b sins) and L is 
the diameter of the major axis of E(s, L), the so called major diameter. 
This is why we require L to be greater than the distance between the 
foci, 2a. Note that E(0, L) = E(2ir, L) and so [0, 2tt] can be viewed as 
the unit circle. 
Let 

Y = {(s,L)\s £ [0,2tt], L > 2a}, 
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then Y is the set of parameters for the ellipses. 

Definition 2.1. Let (s,L) G Y. The elliptical Radon transform of a 
locally integrable function / : R 2 — > R is defined as 

TZf(s,L)= J f(x)dl(x) 

x£E(s,L) 

where dl is the arc length measure on the ellipse E(s, L). 
The backprojection transform is defined for g G C C (Y) as 

K*g{x) = \ g(s, \x - j R (s)\ + \x - j T (s)\)ds 

Js£[0,2it] 

where the positive smooth measure ds is chosen so that 7Z* is the L 2 
adjoint of 1Z. Note that 1Z* integrates over a compact set, [0, 2ir] (or 
the circle), and so 7Z* can be composed with 7Z. Using the parametriza- 
tion of ellipses (s,L) one sees that lZ*g{x) integrates (with a smooth 
measure) over the set of all ellipses passing through x. 

We now briefly discuss the general framework of Guillemin and Stern- 
berg into which our elliptical Radon transform falls. We use this to 
understand the microlocal analysis of 1Z. We begin with some general 
notation we will use when we discuss microlocal analysis. 

2.2. Microlocal Definitions. Let X and Y be smooth manifolds and 
let 

C C T*(Y) x T*{X), 

then we let 

C = {(y,r),x,£)\(y,r),x,-£) G C}. 
The transpose relation is C* C T*{X) x T*(Y): 

C l = {(x,£,y,rj)\(y,ri,x,€) G C} 
If V C T*(X) x T*(Y), then the composition V o C is defined 
VoC = {(x',£,x,OHy,v)eT*(Y) 

with (x',^',y,rj) G V, (y,rj,x,t,) G C}. 

2.3. The Radon Transform and Double Fibrations. Guillemin 

first put the Radon transform into a microlocal framework, and we 
will use this approach to prove Theorem 3.1. These results were first 
given in the technical report [4] (some of which appeared in [5]), then 
outlined in [6, pp. 336-337, 364-365]. The dependence on the measures 
and details of the proofs for the case of equal dimensions were given 
in [13]. Guillemin used the ideas of push- forwards and pullbacks to 
define Radon transforms in [4] and he used these ideas to define Fourier 
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integral operators (FIOs) [4, 6]. Finally Guillemin summarized this 
material in [5]. 

Given smooth connected manifolds X and Y of the same dimension, 
let Z C Y x X be a smooth connected submanifold of codimension 
k < dim(X). We assume that the natural projections 

Z 

"V V* (2-2) 
Y X 

are both fiber maps. In this case, we call (2.2) a double fibration. 

Following Guillemin and Sternberg, we assume that ttr is a proper 
map; that is, the fibers of tir ■ : Z — > X are compact. The double 
fibration allows us to define a Radon transform as follows. For each 
y G Y let 

E{y) = ir R (^({y})), 

then the sets E(y) are all diffeomorphic to the fiber of the fibration 
ttl : Z — > Y. For each x G X let 

F(x) = n L {^({x})) , 

and the sets F(x) are all diffeomorphic to the fiber of ttr : Z — > X. 
Since 7Tr is proper, the sets F(x) are compact. By choosing smooth 
nowhere zero measures on Z, X, and on Y, one can then define a 
smooth nowhere zero measure \x y on E(y) and \i x on and a Radon 
transform 

= / f(x)diMy(x) 

Jx£E(y) 

and the dual transform is 

JyeF(x) 

[6] (see also [13, p. 333]). Since the sets F(x) are compact, one can 
compose R* and R for / G C c pT). If 7r^ is not a proper map, then 
one needs cutoff functions to compose R* with R. We assume E(yi) = 
E(y2) if and only if y\ = y-i and similarly for F . 

Guillemin showed ([4, 5] and with Sternberg [6]) that R is a Fourier 
integral distribution associated with integration over Z and canonical 
relation C = (N*(Z)\{0})'. To understand the properties of R*R, one 
must investigate the mapping properties of C. Let 11^ : C — > T*(Y) 
and II# : C — > T*(X) be the natural projections. Then we have the 
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following diagram: 

C 

n V V* (2-3) 
T*(Y) T*(X) 

This diagram is the microlocal version of (2.2). 

Definition 2.2 ([4, 5]). Let X and Y be manifolds with dim(y) = 
dim(X) and let C C (T*(Y) x T*(X))\{0} be a canonical relation. 
Then, C satisfies the Bolker Assumption if 

ITy : C -> T*(Y) 

is an injective immersion. 

This definition was originally proposed by Guillemin [4], [5, p. 152], 
[6, p. 364-365] because a similar assumption for the Radon transform 
on finite sets implies R*R is injective in this case. Guillemin proved 
that if the measures that define the Radon transform are smooth and 
nowhere zero, and if the Bolker Assumption holds (and R is defined 
by a double fibration for which ttr is proper), then R*R is an elliptic 
pseudodifferential operator. 

In the definition in [4, 13], the dimensions of X and Y are equal, 
but in [5], dim(Y) > dim(X). We use the former definition since, in 
our case, dim(F) = dim(X). Since we assume dim(F) = dim(X), if 
LTy : C — > T*(Y) is an injective immersion, then Uy maps to T*(F)\{0} 
and Tlx is also an immersion [7]. Therefore, Ux maps to T*(X)\{0}. 
So, under the Bolker Assumption, C C (T*(F)\{0}) x (T*(X)\{0}) 
and so R is a Fourier integral operator according to the definition in 
[14]. 

For our transform TZ, the incidence relation is 

Z = {(s,L,x) C Y x D b \x E E(s,L)}. (2.4) 
The double fibration is 

Z 

"V V* (2-5) 
Y B b 

and both projections are fiber maps. Note that the fibers of ttl are 
ellipses, E(s,L) = 7r^(7r£ 1 ({(s, L)}). Furthermore, the fibers of 7Tr are 
diffeomorphic to [0, 2tt] or the circle, so 7Tr is proper. Therefore TZ 
and 1Z* satisfy the conditions outlined above so that Guillemin and 
Sternberg's framework can be applied. 
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3. The Main Result 



We now state the main result of this article. 

Theorem 3.1. Let a > be a constant and let 7t(s) = (cos(s — 
a), sin(s — a)) and 'Jr(s) = (cos(s + a), sin(s + a)) for s £ [0, 2ti] be the 
trajectories of the emitter and receiver respectively. Denote by £'(D&) 
the space of compactly supported distributions supported in the disc, 
of radius b centered at 0, where b = cos a. The elliptical Radon trans- 
form 1Z when restricted to the domain £'{Db) satisfies the microlocal 
conditions in Section 2.3, and 1Z is an elliptic Fourier integral operator 
(FIO) of order -1/2. LetC C T*(Y) x T*(D b ) be the canonical relation 
associated to 7Z. Then, C satisfies the Bolker Assumption (Definition 



As a consequence of this result, we have the following corollary. 

Corollary 3.1. The composition ofTZ with its L? adjoint TZ* when re- 
stricted as a transformation from S'(Db) to ©'(Dj,) is an elliptic pseudo- 
differential operator of order — 1 . 

Remark 3.2. This corollary shows that, for supp / C Df,, the singulari- 
ties of TZ*TZf are at the same locations and co-directions as the singu- 
larities of /, that is, the wavefront sets are the same. In other words, 
1Z*1Z reconstructs all the singularities of /. If D is an elliptic differential 
operator on Y, then one can create an elliptic local reconstruction op- 
erator 1Z*D1Z that will emphasize singularities. REU research student 
Howard Levinson [8] refined and implemented an algorithm of Prof. 
Quinto's of this form and the algorithm does show all the singularities 
of /. Because of the derivative, the algorithm highlights boundaries, 
and this type of algorithm generalizes Lambda Tomography [3, 2]. 

4. Proofs of Theorem 3.1 and Corollary 3.1 

Proof of Theorem 3.1. First, we will calculate C = (N* (Z)\{0})' where 
Z is given by (2.4). Then most of the proof is devoted to showing that 
C satisfies the Bolker Assumption. 
Since Z is defined by the equation 



L-\x- 7r(s)| -\x- 7 T (s)| = 0, 
the differential of the function is a basis for N*(Z). That differential is 



2.2). 
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Note that we are using the fact that 'Jr(s) and 7#(s) are orthogonal 
(as are 7 T (s) and i T (s)). Therefore, C = (N*(Z)\{0})' is 

t * - Jr(s) + x - Ms) \ ^ 
\\x-7 R (s)\ \x-^ t {s)\J 

\(s,L,x) <EZ, u6l\{0} 

The Schwartz kernel of TZ is integration on Z (e.g., [13, Proposition 
1.1]) and so is a Fourier integral distribution associated to C [5]. 

Now we show that the projection 11^ satisfies the Bolker Assumption. 
If A G C the projection is given by (4.1) 

n l( A) = (., L,J-( ^Mf) + As + AL 

V. V. VF-7fl(s)l \x-jr{8)\J 

From this, we have determined, s, L, and u and we need to find x 
knowing that x G E(s, L). 

The easiest way to do this is to develop coordinates, first on the 
ellipses, then on Z, and finally on C. 

We will reduce to the case s = tt/2, so first we give coordinates on the 
ellipse E(ir/2,L) (see Figure 3), which has foci (±a, 6): for (p G [0, 2tt] 
and L > 2a we define 

x' = x'(L, <fi) = (a cosh p cos <p,b + a sinh p sin <fi) 

where p = arccos(2L/a). ^ 

Note that the ellipse E(s, L) meets the ball D& if and only if L G 
(2a, 2\/a 2 + 46 2 ), and for each L G (2a, 2y 'a 2 + 46 2 ), there is an inter- 
val (0 mirl (£),0max(£)) with 71 < 4> mia (L) < max (£) < 2lX SUch that 

x'(L, 4>) G D 6 if and only if G (0 min (L), max (L)). 

Next we let be the counterclockwise rotation about the origin 
through 9 radians. For (s, L) G Y, this gives coordinates on the ellipse 
E(s,L) 

x(s, L, (p) = T s _ 7r / 2 (x'(L, (p)) 

since T s _ 7r / 2 rotates 7^(71-/2) to Jr(s) and rotates 7t(?t/2) to 7r(s). 
This rotation preserves distances and so it does not change the major 
diameter of the ellipse, L. Furthermore, for G (0 min (L), <j> max (L)), 
x(s,L,<t>) G D 6 . 
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Figure 3. Elliptic Coordinate System. The radius of 
the circle in the figure is b. 

We get coordinates on Z as follows. Let 



Z = {(s, L,<f>)\se [0, 2tt], L e (2a, 2Va 2 + 4o 2 ), 

(L),0 max (L))}. 
This provides local coordinates above on Z: 

Z 3 (s, L, (ft) z(s, L, (ft) = (s, L, x(s, L, </>)) . 
Finally, this gives coordinates on C: 

x ■ i R (s) x ■ i T (s) 



A(s, L, 0, lo) = (s, L, u 



2, u 



x-Jr(s)\ \x-j t (s] 
x-Jr(s) x-jt(s) 



+ 



ds + dL ) , 
dx) 



|x-7k(s)| |x-7t(s) 
where (s,L,^)g2x(R\{0}), 
and x = x(s, L, 0) 
and this gives 11^ in coordinates as 

| X " 7 ^l l + . X " 7 ^. 



|x-7 R (s)| |x-7t(s 
Using a series of estimates, we will show the following claim. 



+ ujdLj . 
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Claim 4.1. For (s,L, 0) G Z, the derivative of 



x ( s , L, 0) ■ 7b(s) | x{s,L,<l>)-i T {s) ^ 



\x(s,L,<j>) ~1r(s)\ \x(s,L,<j>) -Jt(s)\ 
with respect to is never zero. 

This claim shows that ITl is injective for the following reasons. Since 
s,L and u ^ are given from the projection 11^, to show IT^ is in- 
jective, all we need to show is that for fixed (s, L) with s G [0, 2n] 
and L G (2a, 2\ / a 2 + 46 2 ), (4.3) determines G (0 min (L), max (L)) 
uniquely. However, Claim 4.1 shows for fixed (s,L) and in the inter- 
val (4>min{L), (f> max _(L)) that (4.3) is either strictly increasing or strictly 
decreasing. Since (0 m i n (L), max (£)) is an interval, this will show that 
(4.3) is an injective function of and therefore IT^ is injective. 

Next, we will use Claim 4.1 to argue that Hl is an immersion. Since 
s, L and oj ^ are given directly from the projection Ul, writing out 
the derivative matrix of Hl shows that one only needs to prove for 
(s, L, 0) G Z and u ^ that the derivative of (4.3) with respect to 
is nowhere zero. This again follows from the claim. 

Proof of Claim 4-1- We prove this by making a reduction to s = tt /2 
and then by a series of estimates. 

Using rotation invariance: we apply T^iis on M 2 , and use the facts 
that distances and dot products are preserved (so e.g., x ■ 7^.(s) = 
x' ■ 7^(7r/2) and \x — Jr(s)\ = \x' — 7 R (7r/2)|) where 

T n/2 - s {x{s, L, 0)) = X'(L, 0), 

Tn/2- s (lR(s)) = 1r{tx/2) = (-a, 6), 
T n/2 _ s (j' R (s)) = 7fl( 7r / 2 ) = (-6, -«)) 
T^/2-si.lTi.s)) = 7t(tt/2) = (a, 6), 

T„ /2 _,(y T (s))=Y T (7r/2) = (-&,a) 
This gives the simplified expression that is equal to (4.3): 

s(l <!>)=( x ' (kL,< ^ ) ' (~ b > + x '( Li ^ ' (~ 6 ' a ^ 



|x'(L,0)-(-a,6)| |x'(L,0)-(a,6)| 
We have reduced showing Ul is an injective immersion to showing 

for L 6 (2o,2V^+46=), e («W,(£),«W£))- 
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Using (4.2) we obtain 



\x'(L, 0) — (—a, 6)| = J a 2 (cosh p cos + l) 2 + a 2 sinh 2 p sin 2 i 
= a (cosh p + cos0), 



\x'{L, 0) — (a, 6)| = J a 2 (cosh pcos — l) 2 + a 2 sinh 2 p sin 2 . 
= a (cosh p — cos 0) 

and 

—6 cosh p cos cf) — b — a sinh p sin 



S{L,<t>) = 
+ 



cosh p + cos 
-6 cosh p cos + 6 + a sinh p sin < 



cosh p — cos 

Adding these two fractions, and simplifying the expression yields the 
following: 

2 sinh p cos (a sin — 6 sinh p) 

S(L, 0) = H{p, 0) = — — — 

(cosh p + cos 0) (cosh p — cos 0) 

Our problem therefore is reduced to showing that J^H(p, 0) ^ 0. 

Denote t = cos0, then sin0 = — t 2 , since G [n, 2tt]. 
Consider a new function H(t) defined as follows (the term —2 sinh p 
does not matter): 

jjf£\ b s ^ nn Pt + atyl — t 2 
cosh 2 p — t 2 

H'(t) = 2 \ 2t ( bt sinh P + at ^ 1 - f2 ) 

(cosh 2 p-i 2 ) L v 1 

+ ^6 sinh p + a\/l - t 2 - (cosh 2 p - t 2 ) . 

Denoting the numerator of the above expression by N, and distributing 
the product we get 



iV = cosh 2 p f 6 sinh p + ay/l — t 2 ; ) 

+ t 2 (b sinh p + ay/1 - t 2 + at ) 



at 2 



(cosh 2 p + t 2 ) (b sinh p + aVl - t 2 ) - (cosh 2 p - t 2 ) 

"(4.4) 
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Denoting 

k = sinh p, we have cosh 2 p = 1 + k 2 , 

and 

N= Vl - a 2 k (1 + k 2 ) + a(l + k 2 )Vl - t 2 + Vl - a 2 k t 2 

o , ^ at 4 a(l + fc 2 H 2 

+ aiVl - t 2 + 



Vi^I 2 VT^¥ 

Notice that all terms in the above expression are non-negative except 
the last one. Hence, in order to show that N > it is enough to show 
that JVi > 0, where 

at 4 a(l + k 2 )t 2 



N x = a(l + k 2 )Vl - t 2 + atVl - t 2 



Since multiplication by \/\ — t 2 does not change the algebraic sign, 
N > if and only if \/l — t 2 N\ > 0, which is equivalent to 

Definition 4.1. We will call the pair (p,4>) admissible if the point 
defined by elliptic coordinates (p, (p) is located inside D(0,b), the disc 
of radius b centered at the origin. 

We need to show that N > for all admissible pairs (p, 4>). 

Lemma 4.2. If \ cos <p\ > b then (p, <fi) is not admissible for any p > 0. 

Proof. Recall the coordinate system (p, <fi) (see Fig. 3) for x = (21,22) 
defined by 

21 = a cosh p cos (AR\ 
%2 = b + a sinh p sin <p ' 

where a > is fixed, p > 0, and <ft £ [0, 27r]. It is easy to notice that 
the coordinate curves (hyperbolas) corresponding to fixed values of 
intersect the disc D(0, b) only for limited values of <fi. To find the range 
of these values consider the following system of equations: 

x\ \ 2 f X2 — b\ 2 kir 



+ — , fc = 0,...,4 



a cos (J) J \a sin J 2 

2 2 +2 2 = b 2 . 

Simplifying the system one gets a quadratic equation with respect to 

(2 2 - 6) 

(22 - b) 2 + 26 sin 2 0(2 2 - b) + a 2 sin 2 cos 2 = 0, 
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which does not have any real roots if | tan0| < ^. Recalling that b = 

\J\ — a 2 we conclude that the hyperbolas corresponding to coordinate 
curves with constant <fi intersect the domain of the function support 
D(0, b) only when \t\ < b = y/1 - a 2 . □ 



4 -3a 2 



Lemma 4.3. //cosh p > — then (p, 0) is not admissible for any 

a 2 



Proof. From equations (4.6) it is easy to notice that the coordinate 
curves (ellipses) corresponding to fixed values of p do not intersect 
the disc D(0, —b,b) if p is large. Namely, the largest value of p for 
which the corresponding ellipse intersects the disc satisfies the equation 
2b = a sinh p. Since cosh 2 p = 1 + sinh 2 p, and b 2 = 1 — a 2 we get that 

l2 4 - 3a 2 

cosh p < — 

a 2 

for each coordinate ellipse that intersects the disc. □ 

Proposition 4.1. If a > 0.8, then for all admissible pairs (p,4>) we 
have N > 0. 

Proof. Recall from Lemma 4.3 that for admissible pairs k 2 + 1 < 
— , hence < k < . Consider the right hand side of 

2 n 



a 

equation (4.5) 

/(*) 



1 + 2P' 

The function f(k) is monotonically decreasing and for admissible pairs 
reaches its minimum 



f(k ■)- k - 2VT ^" 

8 — 7a 2 a 

Hence, if we show that under the conditions of the theorem 

9 4 - 3a 2 , „ „, 

4 < rw (4 - 7) 

then by equation (4.5) and the preceding argument N > 0. At the 
same time by Lemma 4.2 we know that t 2 < 1 — a 2 . So if a > satisfies 

2 4 -3a 2 . 

'- a< rw' (48) 

then (4.7) is satisfied. To finish the proof, notice that the hypothesis 
of the proposition a > 0.8 implies (4.8). □ 
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Proposition 4.2. If a < 1.2, then for all admissible pairs (p, 0) we 
/iowe iV > 0. 

Proof. From Lemma 4.2, we need to consider only those pairs (p, (f)) for 
which |t| < 6, where t = cos0. We will show that iV is positive for this 
range of t. 

The term N in (4.4) can be rewritten as 

N(t) = bVl-t 2 sinh p(cosh 2 p + t 2 ) + a(cosh 2 p - t 2 (cosh 2 p + sinh 2 p)) 

(4.9) 

Using the fact that b > \t\, we have a/1 — t 2 > Vl — = a and we 
substitute this into (4.9) and use the fact that b = yl — a 2 to get 

N{t) > ba sinh p(cosh 2 p + t 2 ) + a(cosh 2 p — t 2 (cosh 2 p + sinh 2 p). 

Factoring a, we are left with wanting to show that 

b sinh p(cosh 2 p + t 2 ) + cosh 2 p — t 2 (cosh 2 p + sinh 2 p) > 0. 

Isolating t 2 , we need to show that 



(fe sinh p + 1) cosh 2 p ^ ^ 
cosh p + sinh p — b sinh p 

(Note that the term cosh 2 p + sinh 2 p — fesinh p is always positive.) We 
show that for b > the left hand side of (4.10) is greater than 1. 
Solving the inequality by let x = sinh p we get the quadratic inequality 

,2 ^ „ • , 1± VI - 86 2 

ox — x + 26 > with roots x = . 

26 

There are no roots if 1 < 86 2 or 1/(2V2) < 6. Since 6 = cos a, this 
corresponds to a < 1.21 radians. □ 

We have now proved Claim 4.1, and from the discussion following 
this claim, we have that LT^ is an injective immersion: the Bolker 
Assumption holds. □ 

Now that we know the Bolker Assumption holds, as mentioned after 
Definition 2.2, the projections 11^ and 11^ map away from 0. Therefore, 
1Z is a Fourier integral operator [14]. Since the measure of integration 
on the ellipses, arc length, is nowhere zero in D{,, 1Z is elliptic. The 
order of 1Z is given by (dim(F) — dim(Z))/2 (see e.g., [5, Theorem 1] 
which gives the order of 1Z*1Z). In our case, Z has dimension 3 and Y 
has dimension 2, hence TZ has order —1/2. 

□ 
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Proof of Corollary 3.1. The proof follows from Guillemin's result [5, 
Theorem 1] as a consequence of Theorem 3.1. We will outline his proof 
that 1Z*1Z is an elliptic pseudodifferential operator since the proof for 
our transform is simple and instructive. 

By Theorem 3.1, TZ is an elliptic Fourier integral operator associated 
with C. By the standard calculus of FIO, TZ* is an elliptic FIO associ- 
ated to C*. Note that we can compose TZ* and TZ because TZ* integrates 
over [0, 27r] (in general, because n R is proper). 

Because the Bolker Assumption holds above C is a local canonical 
graph and so the composition 71*71 is a FIO for functions supported in 
Dfe. Now, because of the injectivity of ITy, C l o C C A where A is the 
diagonal in (T*(Df,)\{0}) 2 by the clean composition of Fourier integral 
operators [1]. 

To show C = A, we need to show ITr : C —> T*(Db)\{0} is surjective. 
This will follow from (4.1) and a geometric argument. Let {x,£) G 
T*(Dfe)\{0}. We now prove there is a (s,L) G Y such that the ellipse 
E(s,L) is conormal to (x,£). First note that any ellipse E(s,L) that 
contains x must have L = \x — 7r(s)| + \x — 7r(s)|. As s ranges from 
to 2-7T the normal line at x to the ellipse E(s, \x — 7r(s)| + \x — 7t(s)|) 
at s rotates completely around 2n radians and therefore for some value 
of So G [0, 2tt] E(sq, \x — 7i?(so)| + \x — 7t(so)|) must be conormal (x,£). 
Since the ellipse is given by the equation L = \x — Jr(s) \ + \x — 7r(s)|, 
its gradient is normal to the ellipse at x; co normals co-parallel this 

gradient are exactly of the form £ = u ( , x ~' yR( } So \ l + a ~ 7T [ s °| ) dx for 

& J ^ \F-7fl(«o)| F-7t(so)|/ 

some to ^ 0. Using (4.1), we see that for this sq, x, oj and L = 
\x — 7_r(so)| + \x — 7t(so)|, there is a A G C with II^(A) = (x, £). 
This finishes the proof that Ur is surjective. Note that one can also 
prove this using the fact that 7Tr is a fibration (and so a submersion) 
and a proper map, but our proof is elementary. This shows that 71*71 
is an elliptic pseudodifferential operator viewed as an operator from 
S'(Db) — > V'(Db). Because the order of TV and 7Z are —1/2, the order 

of n*n is -l. □ 
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